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Abstract 

Recent data in the neutrino sector point towards a relatively large value of the reactor 
angle, incompatible with a vanishing ^13 at about Scr. In order to explain such a result, 
we propose a SUSY model based on a broken U{1)f flavour symmetry with charge 
Le — L^—L-r for lepton doublets and arbitrary charge for the right-handed SU{2) singlet 
fields, where large deviations from the symmetric limit ^12 = vr/4, tan023 ~ 0(1) and 
^13 = mainly come from the charged lepton sector. We show that a description of 
all neutrino data is possible if the charged lepton mass matrix has a special pattern of 
complex matrix elements. 
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1 Introduction 



Hints for not so small 613 in the neutrino sector have been recently confirmed by analyses 
of global neutrino data, including the latest T2K [l] and MINOS [2] results, providing 
sin ^13 > at 3o" CL [3||4]. This has triggered a lot of attention in the community, aimed at 
explaining how such relatively large values can naturally be obtained in a consistent model 
for fermion masses and mixings. Attempts in this direction have been done, for instance, in 
the context of GUT theories [5|, where corrections to the bimaximal or tri-bimaximal mixing 
were responsible for large deviation of ^13 from zero, or using some neutrino mass textures 
giving non- vanishing leading order reactor angle [6]. In this paper we want to analyze the 
possibility to get sizable ^13 from broken U{1)f flavour symmetry [T] with charge Le — L^ — Lr 
for lepton doublets [s] and arbitrary right-handed charges [9]. It is well known that, in the 
limit of exact symmetry, neutrino masses and mixings are described by the following mass 
matrix: 



= mo 1 , (1) 




leading to a spectrum of inverted type, 612 = 7r/4, tan ^23 = x (i.e. large atmospheric mixing 
for X ~ 0{1)) and ^13 = 0. A main problem of such texture is that two eigenvalues have 
the same absolute values, thus preventing the description of the solar mass difference. In 



this context, successful tentatives to solve such problems have been presented in 10 where, 
however, either the solar angle was too large compared to the current best fit point or the 
reactor angle was (almost) vanishing. Corrections of order of the Cabibbo angle A from the 
charged lepton sector [ll|[T3] can simultaneously affect the deviation of 612 from maximal 
mixing and of ^13 from zero, thus allowing a sizable reactor angle, but the question of how to 
naturally explain the smallness of the solar-to-atmospheric mass ratio r (without fine-tuning 
among the model parameters) remained unanswered. In this paper we want to contribute 
to the discussion presenting a model where the previous issues are solved/mitigated. In 
particular, we work in a SUSY framework where the U{1)f flavour symmetry is broken by 
the vevs of two complex fields (j) and 6 of charge = 1 and Q0 = —1/2. We show in Sect|2] 
that an appropriate breaking of the L^. — L^ — symmetry in the neutrino sector is enough 
to guarantee a value of r ~ (and to preserve the leading order (LO) prediction of large 6'23) 
but not to generate the necessary deviations for the solar and reactor angles. The corrections 
from the charged lepton sector are discussed in Sect j3] and their impact on the determination 
of the leptonic mixing angles is investigated in Sect|4| where we show that complex Yukawa 
couplings are needed to make the model in agreement with the experimental data. In Sect. 5 
we draw our conclusions. 



2 The neutrino sector 

The three families of SU{2) lepton doublets U are charged under the non-standard lepton 
number Le — — L^, with values: 

k ^ Le, - L^ - Lr ^ (1, -1, -1) . 

From the dimension 5 operators and in the limit of unbroken Lg — — Lr symmetry, the 
LO lagrangian can be written in the (simplified) form: 

= 77 Melf. + XisUr) H^Hu + h.c. , (2) 
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where M is a large mass scale, Hu represents the (uncharged) Higgs field and Xij are generic 
(9(1) coefficients, that we take as real from now on. The previous lagrangian gives rise to 
the neutrino mass matrix in eq.([T|, where mo = v^xyijlM and x = Xi^/xi2- It is well know 
that such a matrix has two degenerate eigenvalues mi = —m2 = mo y/l + x^ and a vanishing 
one ms = 0; the diagonalizing matrix can be cast in the following form: 



V2 



1 

f 







X 



X 



\ V2Vl+x^ 



1 



\ 



J 



(3) 



12 



from which tan 6 
mass difference Am^^^ 



1, tan 6' 



23 



X and sin6'i3 = 0. From the expressions of the atmospheric 
mip — |m3p, we easily get: 



X 



ml 



-I 



(4) 



so that, for a natural value x ~ 0{1) (needed also to fft the atmospheric angle) we can ffx 
the overall scale mo = 0.035 eV. 

In order to split the solar mass difference and reproduce a correct value of r, we need 
to fill some of the vanishing matrix elements of eq.Q with entries of appropriate order of 
magnitude. To this aim, it would be enough to assume that the U{1)f symmetry is broken 
by the vev of one scalar field; however, as it will be discussed at length in SectjS} we also 
want to generate a large (and avoid large contributions to /i — )■ 67 fro the charged lepton 
sector); this forces to assume that the U{1)f symmetry is broken by the vev of two fields 
(p and 9 of charge = 1 and Qe = —1/2. To allow for non- vanishing vevs, we assume to 
work in the limit of exact SUSY; in this case, the vev of the scalar fields are determined by 
setting to zero the corresponding scalar potential 12 . For the fiavon fields 6 and 0, related 



to the Froggatt-Nielsen symmetry, the non vanishing vevs are determined by the D-term 
associated with the U{1)fn symmetry: Vd = ^{Mpj — — gF\(p\'^)^, where qf is the 

gauge coupling constant of U{1)f and Mpj denotes the contribution of the Fayet-Iliopoulos 
term; the condition for having SUSY minima, Vd = 0, gives the relation for non-vanishing 
vev of the combination Mpj = 5'_f|6'P + It is then natural to assume a common order 

of magnitude for both vevs, ($)/A = A, A being the Cabibbo angle. With this in mind, the 
next-to-leading order (NLO) contributions to the neutrino mass matrix come from: 



C 



NLO 



1 

M 



02 04 \ 



SO that 



m. 




(5) 



(6) 



with Xi = x\jxYi- We see that the charge of is needed to fill the 2-3 sub-block with entries of 
0{}?\ whereas at the same order the contribution of the Q field is negligible. From standard 
perturbation theory, we get the following eigenvalues (in units of mo): 

X2 + x(2a;3 -|- XX4) 



mi 



m2 



ma 



x^ + 



2(l + a;2) 

X2 + X(2X3 -|- XX4 

2(1 +x2) 



A^ 



(7) 



^ X^X2 — 2XX3 -|- X/^ 
1 + X2 



A^ 
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Notice that the corrections to mi^2 are exactly the same. We can compute the small param- 
eter r = Am^g;/Am^^^ in the limit of real parameters, obtaining: 



-2 



2XX3 + X2 + x'^Xi 



'l + x 



2^3/2 



(8) 



For 0{1) parameters, a natural suppression of r by is achieved and the usual problem 
of having r ~ (especially when using non-abelian groups with neutrinos in triplet 

representations) is circumvented here. Notice that the inverted mass ordering predicted 
by the model can be already tested after the second phase of the T2K experiment fT4| or. 



possibly, at future Neutrino Factories or (3 beams 15 



The normalized eigenvectors corresponding to mi and m2 receive 0{\^) corrections 
whereas, at the same order, only the first component of the third eigenvector is modified 
compared to the unbroken limit. It it easy to extract the values of the mixing angles (before 
the diagonalization of the charged leptons): 



sin 6'^3 



tan 



tan 6*23 



x{xX'i + X2- X4) - X3 

{l + x^f^ 



X{XX4 + 2X3) + X2 



2(l + x^ 



,3/2 



1 - 



(9) 



X . 



We see that, although the atmospheric angle is still of (9(1), the solar angle remains too large 
(and always smaller than maximal mixing, due to the condition r > 0) and that the reactor 
angle is as small as A^. The corrections from the charged lepton sector must be large, of 
0{\), to shift both 612 and 613 closer to their recent best fit values [s] while not destroying 
a large 623. Notice that the rotation from such a sector also affects the radiative decays 
of muons and taus. In particular, the correction in the (1 — 2) sector is the most relevant 
one and should be taken under control to avoid a large contribution to // — )■ 67. In fact, it 



has been shown in 16 and 17 that, in a effective theory approach, the amplitudes for the 



radiative decays of the charged leptons depend on the off-diagonal elements of the dipole 
operator Ai, in a basis where the mass matrix of the charged leptons is diagonal. In the case 
of /i — )■ 67, the relevant matrix element is A^i2, suppressed by a cut-off scale and directly 
related to the (1 — 2) charged lepton mixing. In SUSY theories, if we want the cut-off scale 
at the level of ~ 1 TeV, the (1 — 2) mixing should be at least of (9(A^) to avoid a large Aii2 
and then a large fi e'-f branching ratio. This means that the deviations from 612 = vr/4 
can only originate from a different subsector of the charged lepton mass matrix. 



3 The charged lepton sector 



We start with a general consideration about the charged lepton rotation 13 . Let us 



define a generic 3x3 mixing matrix in terms of rotations in three different subspaces: 
/I \ / c?3 sl,e^'\ ( c\2 s\2 0\ 

^ = \^ S\A\ 1 ~S\2 C\2 . (10) 

where we used the superscript e to identify the angles in the charged lepton sector. In the 
left-right basis, the mass matrix is diagonalized by a bi-unitary transformation: 

m, = f/imf[/i. (11) 
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Using for Ul the form in eq.(lO) and disregarding possible CP phases, we obtain: 



'^12'^13^e 
mi = I me(— €235^2 ~ ^12'5l3'523) 
,me( — C^2'^23'5l3 '5l2'523) 



/)>') ( /^'^ I />'>') 

"Vv"-12'-23 ■'12 '='13 "^23/ '-13''23"''t 
"Vv ''23''12''13 ''12'='23/ ''13''23'"'t 



(12) 



The right-handed rotation is not involved in the neutrino mixing; inspired by the minimal 

and by the fact that the CKM matrix is almost diagonal, we 



SU(5) relation mg 

tentatively take Ur ~ 1. Now we want to introduce in appropriate expressions for the 
masses and angles. According to the previous discussion, the needed pattern of corrections 
are guarantee if we take S12 to be smaller than (9(A^) and S13 ~ 0{X) [is]; for the moment, 
we do not specify the value of the atmospheric angle. Finally, we factorize out the r mass 
from rrii and use the mass ratios 

TTT-e : rrin : = : : 1 



to get: 



A^4 


A 


^23-^^ 


■523 


"■523-^ 


„e 
^"23 



(13) 



By construction, this matrix provides a left-handed rotation of type (taking only the leading 
contributions in A): 



Ur 



1 


AS2 


A 




^23 


■^23 


^23-^ 


""^23 


^23 



(14) 



which generates the wanted corrections to the solar and reactor angles while avoiding large 
contributions to — )■ ej. The peculiarity of such a matrix comes from the (2-3) sub-block, 
whose matrix elements are strongly correlated in terms of relative signs and magnitudes: 
in fact, there is only one independent element, dictated by the ^23 angle. If we were able 



to naturally obtain eq.(|l3|), then the neutrino mixing angles, from Upmns = U^U^, would 



have the following expressions (disregarding (9(A^) corrections stemming from neutrino NLO 
corrections): 



^23 "T *23-^ 



A 



'23 + '^23^ 



sin 6'i3 
tan 6*12 
tan 6'23 

•-23 "T '='23"'' 

For 6^23 — )■ (and still x ~ 0(1)), we would get the good results: 

A 



(15) 



'23 



+ cl^x 



sin 9 



13 



72 



tan^i2 ~ 1 + C (y2A) tan023~l^ 



(16) 



However, the matrix in eq.(13) is difficult to implement in models based on U{1)f, which 
only allows to accommodate the order of magnitude of the matrix elements. Moreover, as 

Lt- symmetry, the charges of the 



it will be clear later, in models with broken — L^j_ 
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mass matrix elements {m()23 and {me)^^ are exactly the same and we cannot generate any 
hierarchy among them without invoking additional ingredients. In the following, we then 
look for charge assignments of the right-handed fields which allow to reproduce a structure 



as similar as possible to eq.(13) in terms of powers of A; this is as: 



rrii 




(17) 



We emphasize that a matrix like eq.(17) does not naturally produce mixing angles in agree- 
ment with eq.(15). The operators contributing to the mass matrix are of the form: 



[mi 



M 

A 



A 



:i8) 



where aij and (3ij are positive integers controlling the suppression of each matrix element in 
terms of A, with aij + Pij being the minimum allowed sum (otherwise we would take into 
account higher order combinations of fiavon fields) and a^j are complex 0{1) coefficients. 
is taken to be uncharged. With the U{1)f assignment for the leptons as: 



(1,-1,-1) 

~ {Qe, Qfi, Qt) , 



(19) 



and reminding that 
follows: 



+1 and Qe = —1/2, the charges of the operators in eq.(18) are as 



Q 



nil 



1 + Qe + ail - 

-1 + Qe + a2i 
-1 + Qe + a2i 



/5ii/2 

- /32i/2 

- /32i/2 



1 + ai2 - 

-1 + <5/x + tt22 
-1 + <5/x + ^22 



-/322/2 -I + Qr + a23 - f323/2 
-/322/2 -l + Q. + «23-/323/2 



.(20) 



Notice that, having the fi and r the same L^ — L^j^ — charge, the operators contributing to 
the second and third row of the mass matrix must have exactly the same operator structure 
in terms of fiavon fields or, in other words, the coefficients a and /3 must be equal. 

The choice = 1 is untenable; in fact, in this case (Qmjsa is already invariant (and 
then a23 + ^23 = 0) whereas (Qmjis = 2 + ^13 — /3i3/2 and the minimum choice of charges to 
get {Q 



mi ) 13 



is (ai3, /3i3) = (0, 4); in the mass matrix, this would mean that (m£)i3 would 
be suppressed by A^ compared to {rrii) 33 
One possible solution is to take = 0; in this case (q;23, ^23) 



thus destroying the wanted hierarchy as in eq.(17). 

1,0) and (ai3,/3i3 



(0,2) 



and the last column of eq.(17) is 



(0)/A 
(0)/A 




(21) 



With the same reasoning, we can fix the other charges Qe and Q^ to reproduce eq.(17) and 
then the hierarchy among the charged lepton masses. It is easy to show that the choice 
iQcQfi) = (^7,5/2) gives the correct result: 



m£ ~ A 



(0)VA' 
(0)VA^ 






^A^ 


A6 






A^ 


A2 


i) 






A2 





(22) 



The set of charges that can reproduce the observed charged lepton masses and a structure 
similar to eq.(13) are summarized in Tab.([l|. 
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le 


I, 


Ir 













1 


-1 


-1 


-7 5/2 


1 -1/2 



Table 1: Set of charges that can accommodate the observed charged lepton masses and 
neutrino mixings. 



4 The UpMNS mixing matrix 

In the language of U{1)f flavour models, the mass matrix produced by our charge as- 
signment is as follows: 



OiiiX^ (ii2X^ di^X 
mt^mr \a21X'' a22e^'^22^2 ^23 e*<^23 



(23) 



where we have explicitely shown the phases of the (22), (23) and (32) matrix elements (phases 
into the other entries are unimportant or even irrelevant for our reasoning); all ttjj are generic 
0(1) coefficient^ For our discussion, it is important to quote the leptonic (obtained 
after diagonalizing mim\) and the three neutrino mixing angles, including the corrections 
from the charged lepton sector: 



tan 6* 



12 



022 '^23 + ^132 + 2 ^22 ^23 0-32 COs( 

+ 4^0 



^23 



sin^ 



13 



tan 6' 



12 



ai3 



1 + 



A 



'^22 + '^32 + 2 a22 ^32 X COS 0^ 
(1 +X2)D 

2ai3 [a22(a22a^ - ^32(^233: cos(023 - 0a) + cos 0a)) + a23ai2 cos ^23] 



(24) 



A 



tan 9 



23 



033 + 2a23a; cos 023 + 



1 -|- a^gX^ — 2a23X cos 023 ' 
where we have used the short-hand notations: 



D 



''22 



2022023^32 COs(023 



^a) + a23«32 



^32 



In spite of their apparent complexity, the previous relations offer a simple way to reconcile 
the model with the experimental data. First of all, we observe that 912-, ^13 and the deviation 
of 9 12 from maximal mixing are all of 0{X)\ to achieve a suppression for the leptonic 9 12 
we need cos(023 — 0a) ~ —1; on the other hand, this is not enough to guarantee an almost 
maximal atmospheric mixing, unless 023 and x are both very different from 0{1)\ we then 
require the additional condition cos 0a ~ 0. With these positions, the analytical structure 
of the mixing angles is simple and allows a clear understanding of the outcome of such a 



^We note that the mass matrix mi imphes a right-handed charged lepton rotation very close to the 
identity, as assumed in Sect. 3. 
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procedure: 



tan6'^2 



sin 6i3 
tan 6*12 



~ (^22 (^23 + •^32! 



^22 + 023^32 



^13 



2 

23 



' '^22 "I" '^32 •^'^ 



1 + 



0.22 + Ct23 0-32 J \l 1 + 

2 ai3 022 a; 



(25) 



(022 + 023 O32) a/TT 



A 



tan 6' 



O23 + ^ 



23 



1 + ^ 



Bearing in mind that the parameters 022,23,32 are all positive, we see that a cancellation in 
the numerator of tan6'^2 occurs; moreover, we expect: 



sin 6* 



A 



13 



tan 6 



12 



1 ± 



A 



tan 6 



23 



2 ^ " 

and then a clear correlation among the solar and reactor angles: 



e 



12 



~ — ± 

4 



(26) 



(27) 



It is interesting to observe that, in the limit of real Yukawa parameters, it would be very 
difficult to simultaneously take under control both 6^2 and 6'23, unless strong cancellations 
among the relevant Ojj are invoked. Our analytical results are confirmed by a numerical 
evaluation of the mixing angles, obtained extracting randomly the values of the parameters 
Ojj and x,Xi of the mass matrices (23) and ([6]), respectively, in the interval [1/3,3] and the 
phases (023, 0^) around (37r/2, 7r/2) ±0.1; we do not impose any other constraints. In Fig.([T]) 
we show the distributions of values for tan ^^2! we see that many points fall in the region 
below A^ (indicated with a solid vertical line), corresponding to almost 35% of the total 
number of realizations. 




0.00 



0.00 0.05 0.10 0.15 0.20 0.25 
Tan 9U 



0.30 



Figure 1: Distribution for tan (arbitrary units). The black line indicates the value A^ 



In Fig. ([2]) we present the main result of this paper, namely the correlation among 612 and ^13 
in the (sin^ ^13, tan 6*12) plan^ In the central plot we show the simultaneous determinations 

•^We only show values of tan 6*12 < 1; an almost symmetric region is populated at tan0i2 > 1. 
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of both angles; vertical solid (black) lines enclose the 3a bounds on ^^13 from |4] whereas the 
gray horizontal bands are the regions excluded by the experimental data on 612 (also at 3cr 
from 111), that is 



0.001 < sin^^ia < 0.035 



0.61 < tan ^12 <0.75. 



(28) 



We clearly see that the correlation shown in eq.(27) is in fact realized in our model: large 




Figure 2: Central panel: scatter plot in the variables (sin^ ^13, tan^i2) as obtained from our 
model. The gray bands represents the regions excluded by the experimental data on tan 6^12 at 
3a ^ whereas the vertical lines enclose the 3a range on sin^^i3. See text for further details. 
Upper and left panels: distributions of the variables sin^ ^13 and tan ^12, respectively, and 
their 3a allowed regions, enclosed into the solid (black) lines. 



values of the reactor angle bring the solar one into the experimentally allowed region whereas, 
for small 6*13, 6*12 is still close to maximal mixing. In the upper and left panels we also display 
the distributions of the related variables sin^ ^13 and tan6'i2, respectively, with the vertical 
solid (black) lines enclosing their Sex bounds. It is interesting to observe that ~67% of the 
points in the sin^ ^13 distribution falls in the current allowed range whereas such a number 
decreases to ~30% for tan6'i2. We also notice that the distributions of the angles have the 
picks in the region around the values indicated in eq.(26). For the atmospheric angle we 
found no significant correlations with the other mixing parameters and a distribution of 
values with more than 50% in the allowed 3cr range. We do not show the corresponding 
plots. 

Finally, we analyze the model predition for the effective mass appearing in the 
neutrinoless double (3 decay amplitude; since the lightest neutrino mass is at most of 
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(9(A^) (and S13 ~ 0{X)), we can safely neglect its contribution to the effective mass niee 
appearing in the neutrinoless double (3 decay amplitude, which then reads: 

Imeel ~ (4 + sl^ e'n , (29) 

where a is the Majorana phase associated to m2. Since we choose to work with real param- 
eters in the neutrino mass matrix, a is generated by the charged lepton rotation. In a basis 
where they are diagonal, the condition {(1)23, 4>a) = (3vr/2,7r/2) gives a = 0; however, like in 
the case for the Dirac CP phase, when we perturb the previous conditions any value of a in 
[0, 27i] can be generated and we get: 

|mee| ~ (1.9 -5.1) X lO^^el/. (30) 



5 Conclusions 

In this paper we presented a model for lepton masses and mixings based on a U{1)f 
flavour symmetry with charge Lg — L^ — Lr- The symmetry is broken by the vevs of two scalar 
fields and, thanks to the corrections from the charged lepton sector, we get a large 613 ~ 0{X), 
while the other mixing angles can be made compatible with the experimental data. For this 
to happen, we used the freedom given by the complex phases of the Yukawa couplings 
defining the charged lepton mass matrix; we observed that, under certain hypothesis, we 
are able to maintain the atmospheric angle close to maximal mixing and to lower the value 
of the leptonic 12 rotation at least at (9(A^), necessary to avoid dangerous contributions 
to the rare decays of the charged leptons. We also obtained a naturally small value of the 
solar-to-atmospheric mass difference ratio r, at the level of A^. 
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